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Body  of  the  Report 


1 .  Background  of  the  Problem 

It  is  well  known  that  mixtures  of  cyclic  nitramine-polymer  binder  offer  many  advantages  for 
advanced  propulsions.  They  have  better  performances  (large  amount  of  gas,  high  specific  impulse  for  rockets 
and  impetus  for  guns),  safety  (low  sensitivity)  and  environment  friendness  (nontoxic  and  noncorrosive 
combustion  products).  Energetic  binders  are  preferable  for  that  propulsions.  The  study  of  combustion 
mechanism  of  the  nitramine-energetic  polymer  binder  mixtures  is  very  important.  The  understanding  of  the 
mechanisms  is  useful  for  many  applications  of  that  propulsions. 

Different  mixtures  based  on  new  energetic  binders  were  studied  in  this  work.  Especially  detailed 
were  studied  mixtures  with  improved  characteristics. 

2.  Approach 

Different  types  of  energetic  binders  were  used  during  the  contract  period.  They  were: 
glycidyl  azide  polymer  with  molecular  mass  of  two  thousand  (GAP2),  copolymer  of  two 
thermoplastic  elastomers  -  poly  bis  3,3-azidomethyloxetane  (BAMO)  and  poly  3-azidomethyl-3- 
methyloxetane  (AMMO),  in  proportion  50:50,  copolymer  of  BAMO  with  tetrahydrofuran  (THF), 
in  proportion  50:50,  and  glycidyl  azide-polyurethane  copolymer  (GAP1U).  Cyclic  nitramines  HMX 
and  RDX  are  used  as  oxidizers.  The  ratio  20:80  for  components  of  energetic  binder/nitramine  was 
used  for  all  of  the  mixtures.  Nitramine  powders  have  the  following  characteristics:  50%  by  weight 
have  sizes  <50  microns  and  50%  -  150-300  microns.  Purity  of  HMX  and  RDX  was  99%.  The 
propellant  samples  were  0.8cm  in  diameter  and  2  cm  in  length. 

Thus  the  following  mixtures  were  prepared  and  investigated:  1)  GAP2/HMX,  2)  GAP2/RDX, 

3)  (BAMO-AMMO)/HMX,  4)  (B AMO- AMMO)/RDX,  5)  (BAMO-THF)/HMX, 

6)  (B AMO-THF)/RDX  and  7)  GAP1U/HMX,  8)  GAP1U/RDX. 

The  study  of  the  combustion  mechanism  of  the  mixtures  was  performed  by  obtaining  geometrical 
and  thermal  structure  of  burning  waves  of  the  mixtures  and  by  processing  of  the  obtained  data.  The 
peculiarity  of  the  work  is  the  measurement  of  all  characteristics  of  combustion  waves  under  change 
of  external  parameters  -  pressure  and  initial  sample  temperature.  This  approach  allows  one  to 
established  distinctive  features  of  the  combustion  mechanism  of  the  studied  mixtures. 

The  most  reliable  data  were  received  for  mixtures  1)  GAP2/HMX,  2)  GAP2/RDX  and 
3)  (BAMO-AMMO)ZHMX,  4)  (BAMO-AMMO)/RDX.  Because  of  that  temperature  and  pressure 
sensitivities  of  burning  rate  and  surface  temperature  were  obtained  for  these  mixtures.  Burning  rate 
response  functions  to  oscillatory  pressure  were  estimated  by  using  these  sensitivities. 

3.  Experimental  Techniques. 

3. 1.  Conditions  of  experiments  and  microthermocouple  techniques. 

Temperature  profiles  of  the  combustion  waves  and  the  burning  surface  temperatures  were  obtained  by 
microthermocouple  methods.  Profiles  of  the  combustion  waves  were  obtained  by  inicrothermocouples  imbedded  into 
solid.  Thermocouples  went  through  combustion  waves  when  the  waves  propagated  through  the  solid  samples  and 
registered  temperature  profiles.  The  ribbon  U-shaped  thermocouples  made  of  alloys  W+5%Re/W+20%Re  of  2-7mkm 
thick  were  imbedded  into  the  samples.  Cylindrical  samples  of  8  mm  in  diameter  and  3  cm  in  length  were  cat  by 
small  thin  knife  of  U-shape  form,  thermocouple  placed  into  the  created  cross-sections,  which  were  moistend  by 
acetone,  and  the  joined  sections  with  thermocouples  inside  were  protected  by  thin  layer  of  polymer  glue.  Sample 
drying  was  performed  during  6  ours.  Every  sample  had  inside  2-3  thermocouples  placed  one  above  the  other. 
Distances  between  the  junctions  were  2-4  mm.  The  samples  were  burned  in  a  bomb  of  constant  pressure  in 
atmosphere  of  nitrogen  at  pressures  5-100atm  and  at  sample  initial  temperatures  T„=  -100,  +20.  and  +90°C.  In 
experiments  at  elevated  sample  temperatures  the  sample  was  placed  into  a  small  thermostat  inside  the  bomb.  The 
sample  heating  was  controlled  by  thermocouples.  Samples  were  ignited  by  electrically  heated  wire.  Thermocouple 
signals  were  registered  by  amplifier  and  oscillograph.  Burning  rate  was  measured  by  time  delay  between  the 
thermocouple  signals,  by  photoregistrations  of  sample  combustion  and  by  pressure  increase  during  the  sample 
combustion.  Photoregistration  of  sample  combustion  allows  also  the  type  of  combustion  to  be  established.  Burning 
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surface  temperatures  were  measured  by  thermocouples  that  are  being  pressed  to  the  surface  during  sample 
combustion  and  by  establishing  the  locations  of  slope  breaks  on  temperature  profiles  registered  by  thermocouples 
(see  below). 

3.2.  Temperature  profile  measurement  validations.  As  a  rule,  temperature  gradient  close  to  the  burning 
surface  has  a  very  high  value.  It  implies  that  thermocouple  measurements  can  give  temperature  profiles  with 
significant  errors  due  to  thermocouple  heat  inertia.  Because  of  that  it  is  necessary  to  find  conditions  under  which 
thermocouple  measurements  in  combustion  waves  will  introduce  small  errors.  These  conditions  have  been  found  by 
numerical  simulations.  The  thermocouple  partially  absorbs  the  heat  of  the  thermal  layer  and  decreases  the 
temperature  at  the  point  of  the  measurement.  The  requirement  of  small  temperature  error  (less  than  10%)  of  the 
thermocouple  is  indicated  by  the  following  formula:  h<0.2 

Here:  h  -  thermocouple  thickness  (in  cm),  %-  heat  diffusivity  of  the  solid  (in  cm2/s),  -  linear  burning  rate  (in  cm/s). 

X/rb=/  where  /  is  the  thickness  of  heat  layer  of  the  condensed  phase.  There  is  another  requirement  for  correct 
measurements:  the  thermocouples  have  to  have  U-shape  form.  It  is  necessary  because  of  a  very  high  difference 
between  heat  conduction  coefficients  of  metallic  thermocouple  wire  and  that  of  solid  or  gas.  Junctions  of  U-shaped 
thermocouples  do  not  experience  large  temperature  decrease.  Modelling  experiments  and  numerical  simulations 
show  that  the  decrease  of  junction  temperature  will  be  small  (<3%)  if  the  horizontal  part  of  the  U-shaped 
thermocouple  is  about  one  hundred  times  more  than  thermocouple  thickness  h.  Thermal  inertia  of  the  thermocouple 
in  the  gas  phase  can  be  taken  into  consideration  and  eliminated  by  a  correction  procedure.  The  procedure  implies  the 
use  of  the  following  equation: 

dTex/dx=(rbx0)-'-(T  -  Tcx); 

Here:  T  -  the  real  temperature  of  gas  in  the  combustion  wave;  Tex  -  the  temperature  registrated  by  thermocouple;  x„  - 
time  response  of  the  thermocouple  in  gas.  Time  response  depends  on  mass  burning  rate  m  (m=p-rt,)  and 
temperature  T.  The  temperature  profiles  in  gas  were  corrected  by  this  equation.  The  theory  of  thermocouple 
measurements  in  combustion  waves  of  solids  was  created  and  confirmed  by  measurements  of  combustion  wave 
temperature  profiles  by  thermocouples  with  sequentially  decreased  thickness  (method  of  “zero  diameter”).  All  the 
above  mentioned  requirements  have  been  met  in  the  investigations.  Different  types  of  metal  wires  for  thermocouples 
were  used  (We,Re  and  Pt,Rh)  to  test  the  catalytic  effect  on  thermocouples.  The  catalytic  effect  was  not  observed.  The 
method  of  burning  surface  temperature  measurement  by  determining  locations  of  slope  breaks  on  temperature 
profiles  registered  by  thermocouples  (method  of  “slope  break”)  is  based  on  the  existence  of  the  delay  (rb  x„)  on  the 
temperature  profiles  when  thermocouples  go  through  the  burning  surface.  The  delay  is  due  to  the  change  of  heat 
exchange  between  environment  and  thermocouple:  contact  heat  exchange  in  solid  is  replaced  by  convective  heat 
exchange  in  gas. 

4.  Results  of  Measurements  and  Processing  for  Mixtures  1  -4 

As  it  was  mentioned  above,  four  new  mixtures  -  1  -  4  were  detailed  studied.  The  mixtures 
had  more  uniform  contents  which  allow  decreased  scattering  of  all  characteristics  to  be  obtained. 
The  samples  had  increased  solid  density  p  and  improved  mechanical  properties.  The  thermodynamic 
calculated  adiabatic  flame  temperatures  Ta(i  are  equal  to  2595-2787  K. 

The  characteristics  of  the  initial  mixtures  were  as  follows: 

1)  GAP2/HMX,  p=l,52  g/cm3;  sum  formula:  O24.156H21.8i1N26.886C17.896;  Tad=2776  K; 

2)  GAP2/RDX,  p=l,51  g/cm3;  sum  formula:  O24.156H21.8i1N26.886Cn.896;  Tad=2787  K; 

3)  [BAMO-AMMOJ/HMX,  p=l,  63g/cm3;  sum  formula:  O23.88H32.24N26.37C1s.06;  Ta(|=2595  K; 

4)  [BAMO-AMMO]/RDX,  p=  1,56  g/cm3;  sum  formula:  O23.88H32.24N26.31C1g.06;  Tad=2606  K; 

Carbon  residue  on  the  burning  surface  of  the  mixtures  did  not  observed. 

4.1.  Temperature  Profiles  and  Burning  Wave  Characteristics. 

Temperature  profiles  for  the  mixtures  were  obtained  at  pressures  5-100  atm  since  at  p<5  atm  the 
mixtures  can  not  burn  stable  and  their  ignition  is  very  difficult.  The  obtained  temperature  profiles 
have  stable  character,  show  a  small  date  scattering,  and,  as  a  rule,  have  not  regular  pulsations  (only 
at  5  atm  some  of  samples  have  irregular  temperature  pulsations  of  ~  200°  in  amplitude).  10-12 
temperature  curves  were  received  for  every  regime  of  combustion  and  then  averaging  of  these 
curves  was  performed.  Thus  every  regime  of  combustion  presented  below  is  outlined  by  one 
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averaged  temperature  profile.  Figures  1,  2  show  the  obtained  averaged  temperature  profiles  of 
mixtures  1,  2  (GAP2/HMX  and  GAP2/RDX)  and  mixtures  3,  4  [(BAMO-AMMO)/HMX  and 
(BAMO-AMMO)/RDX]  at  pressure  5  atm  and  sample  initial  temperature  To=20°C.  Figures  3-6 
show  the  obtained  averaged  temperature  profiles  of  mixture  1  at  pressures  p=10,  20  50  and  100 
atm  and  at  To=-100,  20  and  +90°C.  Figures  7-10  show  the  averaged  profiles  for  mixture  2,  Figs. 
11-14  -  for  mixture  3,  and  Figs.  15-18  -  for  mixture  4,  under  the  same  conditions.  The  presented 
Figs,  show  that  the  gas  phase  of  the  combustion  waves  of  the  investigated  mixtures  is  comprised  of 
one  or  two  zones.  The  first  zone  of  the  two-zone  structure  is  a  relatively  low-temperature  zone 
with  a  zone  temperature  Ti,  and  the  second  zone  is  the  flame  zone  with  the  final  temperature  Tf.  At 
elevated  pressures  (at  p>20  atm)  the  two  zones  always  merged  into  one  zone.  Adiabatic  flame 
temperature  Tad  is  achieved  at  about  50-80  atm  under  the  investigated  conditions. 

Tables  1,  2,  3  and  4  show  the  averaged  burning  wave  parameters  obtained  by  experiments  and  by 
processing  experimental  data  for  every  mixture  (see  below).  Small  corrections  of  the  experimental 
data  were  performed  during  preparations  of  the  Final  Report. 

4.1.1.  Mass  burning  rates.  Tables  1-4  show  that  mass  burning  rates  m  of  all  mixtures  increase 
when  pressure  p  and  sample  initial  temperature  T0  increase.  Values  of  m  for  mixtures  1-4  are  close 
each  other,  but  the  mixtures  with  RDX  have  a  bit  more  values  of  m.  Mass  burning  rates  of  the 
investigated  mixtures  increase  from  0.10-0.12  g/cm2s  at  5  atm  up  to  1.2  1.5  g/cm2s  at  100  atm. 
Values  of  m  of  the  investigated  mixtures  at  10  atm  increase  from  0.21-0.22  g/cm2s  at  To=20°C  up 
to  0.2-0.28  g/cm2s  at  To=+90°C  and  decrease  to  0.15-0.16  g/cm2s  at  To=-100°C.  Difference 
between  values  of  m  at  different  T0  increases  with  pressure.  For  instance,  at  80  atm  m=0.94-l.l 
g/cm2s  at  To=-100°C,  and  m=l. 25-1.42  g/cm2s  at  To=+90°C.  Standard  deviation  of  mass  burning 
rate  measurements  is  &m=±5%. 

4.1.2.  Burning  surface  temperatures.  Tables  1-4  show  that  burning  surface  temperature  Ts 
always  increases  when  T0  and  pressure  increase.  Values  of  Ts  for  mixtures  1-4  are  also  close  each 
other  and  increase  from  330-335°C  at  5  atm  up  to  492-500°C  (To=20°C)  at  100  atm.  Difference 
between  values  of  Ts  at  different  T0  decreases  when  pressure  increases.  For  instance,  at  10  atm  Ts  is 
equal  to  345-350°C  at  To=-100°C  and  to  375-380°C  at  To=+90°C,  but  at  80  atm  Ts  comprises  475- 
480°C  at  To=-100°C  and  490-500°C  at  T„=+90°C  Standard  deviation  of  burning  surface 
temperature  measurements  is  QTS=±5%. 

4.1.3.  Geometric  structures  of  temperature  profiles  in  gas.  As  it  was  mentioned  above,  the  gas 
phase  of  the  burning  waves  of  the  studied  mixtures  has  two  type  of  the  structure:  single-zone  and 
two-zone  structures.  The  mixtures  always  have  two-zone  structure  at  p<  10-20  atm  and  both  zones 
merge  at  20-50  atm.  Values  of  mean  temperature  of  the  first  flame  Ti  are  equal  to  900-980°C  at  5 
atm  and  increase  up  to  1000-1 100°C  at  20  atm.  Influence  of  T0  on  Ti  is  very  week.  Please,  note 
that  mark  -  in  Tables  1-4  implies  that  the  parameter,  Ti  in  this  case,  does  not  exist  in  this  regime. 
Flame  temperatures  Tf  are  small  at  5  atm  (1800-1950°C)  and  increase  when  pressure  increases  up 
to  thermodynamically  calculated  values  of  the  temperatures  at  50-80  atm.  As  a  rule,  values  of  Tf 
increase  approximately  on  the  value  of  (To-20°)  when  sample  temperature  changes. 

Distances  Li  between  the  burning  surface  and  the  beginning  of  the  flame  zone  decrease  when 
preesure  increases.  In  fact,  Lj  is  the  length  of  the  first  flame.  Values  of  Li  are  equal  to  300-800  pm 
at  5  atm  and  decrease  up  to  100-200  pm  at  20  atm.  Sample  initial  temperature  T0  significantly 
influences  on  values  of  Lp  increasing  of  T0  prolongs  Li  and  can  creates  two-sone  structure  in  the 
case  in  which  at  smaller  T0  was  one-zone  structure.  For  instance,  mixtures  1  and  2  at  To=20°C  have 
not  two  zones  at  p>10  atm,  but  at  T„=+90°C  two-zone  structure  came  to  existence  at  20  atm. 

Tables  and  Figs,  also  show  that  at  10  atm  values  of  Li  of  all  mixtures  increase  1.5-2  times  when  T„ 
increases  from  20°  up  to  +90°C. 
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Distances  L  between  the  burning  surface  and  the  flame  (up  to  0.99  Tf)  also  decrease  when  preesure 
increases  In  fact,  L  is  the  length  of  the  gas  phase  reaction  zone.  Values  of  L  are  equal  to  2-2.5  mm 
at  5atm  and  decrease  up  to  1.0-1. 5  mm  at  100  atm.  As  a  rule,  increasing  T0  increases  L  and 
decreasing  T0  decreases  L.  The  decreasing  can  be  very  significant:  for  instance,  values  L  of  mixture 
4  decrease  about  2  times  at  10  atm  when  T„  decreases  from  20  up  to  -100°C. 

Standard  deviations  of  thickness  L,  measurements  are  ±(20-40)%  and  that  of  L  are  ±(10-20)%. 


Table  1 :  Burning  wave  parameters  of  mixture  1  (GAP2/HMX) 


p,  atm 

5 

10 

20 

50 

80 

100 

No 

oH 

O 

o 

20 

-100/20/90 

-100/20/90 

-100/20/90 

-100/20/90 

-100/20 

1 

m, 

g/cm2s 

0.1 

0.155/0.21/ 

0.28/0.37/ 

0.64/0.77/ 

0.94/1.08/ 

1.1/1.26 

0.26 

0.43 

0.88 

1.25 

2 

O 

O 

330 

345/365/375 

380/400/410 

445/460/470 

475/483/490 

485/492 

3 

(p-10"4,  K/cm 

3.5 

2.0/4.0/5.0 

3.5/5.0/7.8 

9.0/10/11 

11/11.7/12.4 

11.5/12 

4 

q,  cal/g 

64 

35/37/42 

29/31/43 

35/30/33 

29/28/26 

27/25 

5 

Qr>  cal/9 

14 

13/8/9 

8/7/8 

5/5/5 

4/3. 7/3.4 

3.4/33 

6 

Q,  cal/g 

65 

140/112/82 

174/133/96 

212/160/134 

216/173/150 

223/180 

7 

/,  nm 

167 

150/107/106 

108/77/74 

50/58/53 

53/49/49 

48/48 

8 

lm, 

45 

30/40/55 

30/35/40 

40/30/38 

20/30/35 

20/20 

9 

^*103,  cm2/s 

1.2 

1.2/1. 5/1. 8 

2/2/22 

2.7/2. 9/3. 2 

33/3.5/3.8 

3. 5/3. 8 

10 

Xs-io4,' 

7.3 

7.3/9.1/10.8 

12/12/13 

16.4/18/19.4 

19.4/21.2/23 

21.2/23 

11 

cal/cm-s-K 

Ti,°C 

980 

-/1050/1000 

-/-/1 100 

-/-/- 

-/-/- 

-/- 

12 

Tf,  °C 

1900 

1950/2050/ 

2150 

2150/2200/ 

2280 

2400/2500/ 

2550 

2400/2500 

2550 

2400/ 

2550 

13 

Li,  nm 

600 

-/500/900 

-/-/350 

-/-/- 

-/-/- 

7- 

14 

L,  mm 

2.0 

1.5/2.3/2.3 

1.5/1. 6/2.3 

2.0/2.0/2.5 

1.7/1. 5/2 

1.5/1. 7 

15 

9,  jam 

43 

43/24/19 

20/15/13 

21/8.0/7.0 

6.5/53/4.9 

5.6/4. 9 

16 

Xg-104, 

2.0 

2.1/2.2/2.2 

2.3/23/2.4 

2.5/2.6/2.6 

2. 6/2. 6/2.6 

2. 6/2.6 

17 

cal/cms-K 
(Do,  kcal/cm3s 

1.4 

1.3/3.5/5.5 

4.1/7.8/14 

24/32.7/41 

44/53.7/66 

53.8/64 

Please,  note:  mark  -  in  the  Table  implies  that  the  parameter  does  not  exist  in  this  regime. 


4.1.4.  Characteristic  sizes  in  gas  and  solid. 

Characteristic  thickness  /  and  /„,  in  solid  presented  on  Tables  1-4  were  obtained  by  temperature 
distributions  in  solid.  Experimentally  obtained  thickness  of  heat  layer  in  solid  /  is  a  distance  between 
the  burning  surface  and  a  section  in  solid  with  temperature  T*=(Ts-T0)/e  -  T0  (e  is  the  base  of  In). 
Thickness  of  partly  melt  layer  (with  melt  nitramine)  in  solid  lm  is  a  distance  between  the  burning 
surface  and  a  section  in  solid  with  RDX  or  HMX  melt  temperatures  Tm  which  are  equal  to  200  C 
for  RDX  and  280°C  for  HMX.  Values  of  /  and  decrease,  as  a  rule,  when  pressure  increases,  and 
/>/  Conductive  sizes  3  of  the  gas  phase  were  obtained  by  formula:  0=Vmcpi  where  \  is 
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coefficient  of  heat  conduction  in  gas  and  cp  is  coefficient  of  specific  heat  of  the  gas  phase.  Values  Xg 
and  cp  were  estimated  by  content  of  combustion  products.  The  obtained  dependencies  of  ^g(Ts)  and 
Cp(Ts)  are  presented  on  Tables  1-4.  The  Tables  show  that  values  &  decrease  when  pressure 
increases  and  they  are  equal  to  40-51  pm  at  5  atm  and  to  4-5  pm  at  100  atm.  It  is  obvious  that 
&«Li  and  &«L. 


Table  2:  Burning  wave  parameters  of  mixture  2  (GAP2/RDX) 


p,  atm 

5 

10 

20 

50 

80 

100 

No 

T0)  °C 

20 

-100/20/90 

-100/20/90 

-100/20/90 

-100/20/90 

-100/20 

1 

m, 

g/cm2s 

0.11 

0.15/0.21/ 

0.27 

0. 3/0.4/ 

0.47 

0.65/0.84/ 

0.96 

0.98/1.25/ 

1.42 

1.2/1.51 

2 

Ts,  °C 

330 

345/365/380 

385/405/420  . 

450/468/478 

480/492/500 

490/500 

3 

(p-10"4,  K/cm 

2.0 

3.0/3. 5/4. 5 

4.5/5.0/7.0 

8.5/9.5/12 

11.2/12.5/ 

14.7 

12.8/14 

4 

q,  cal/g 

38 

47/37/37 

35/29/36 

31/28/31 

29/26/27 

28/24 

5 

qr,  cal/g 

15 

13/10/9 

7/6. 3/6. 3 

5.3/4.6/4.8 

3. 8/3. 2/3 

3.3/3 

6 

Q,  cal/g 

95 

141/113/93 

176/142/113 

207/170/142 

223/183/157 

228/190 

7 

/,  pm 

165 

120/106/100 

100/85/72 

63/54/50 

52/48/45 

48/40 

8 

lm ,  nm 

110 

55/75/100 

55/65/70 

45/50/68 

40/40/40 

40/38 

9 

X-103,  cm2/s 

1.2 

1.2/1.5/1.8 

1. 8/2/2. 2 

2.7/3.0/3.2 

3. 3/3. 8/4.0 

3.8/4. 2 

10 

Xs-104, 

cal/cm-s-K 

7.2 

7.3/9.0/10.3 

11/12/13 

16/18/19 

20/23/24 

23/25 

11 

Ti,°C 

950 

850/950/ 

1050 

-/-/1100 

77- 

77- 

7- 

12 

Tf)  °C 

1950 

2050/2100/ 

-2180 

2100/2200/ 

2300 

2400/2500/ 

2580 

2420/2515 

2580 

2420/ 

2515 

13 

L[,  nm 

700 

400/500/600 

-/-/200 

77- 

77- 

-/- 

14 

L,  mm 

2.2 

1.5/1. 8/2. 2 

1.7/2.0/2.0 

1.5/2.0/2.5 

1.5/1. 6/1. 5 

1. 2/1.2 

15 

9,  jrm 

46 

34/25/19 

17.5/14/12 

9.0/7.0/6.0 

6.0/5. 0/4.0 

5.0/4. 0 

16 

A-g'104, 

cal/cm-s-K 

2.1 

22/22/22 

2. 2/2. 3/2.4 

2.4/2. 5/2.5 

2.5/2. 5/2. 6 

2. 5/2. 6 

17 

Oo,  kcal/cm3s 

0.9 

2.0/3. 1/5.1 

5.7/8.4/13.8 

23.2/34/49 

47/66/89 

65/90 

4.1.5.  Coefficients  of  heat  dijfusivity  and  heat  conductivity  in  solid  heat  layer. 

The  coefficients  of  solid  heat  diffusivity  were  obtained  by  formula  %=l-r t>-  Tables  1-4  show  that  x 
increases  when  pressure  increases  -  from  1.2-1  O'3  cm2/s  at  5  atm  up  to  4- 1 0‘3  cm2/s  at  100  atm.  That 
are  normal  values  of  x  for  solid  heat  layer  of  many  propellant  mixtures.  Values  of  x  increase  when 
T0  increases.  The  coefficients  of  solid  heat  conduction  were  obtained  by  formula  7.s=x-c-p.  It  was 
assumed  that  solid  specific  heat  C=0.4  cal/g-K.  Tables  show  that  Xs  increases  with  pressure  -  from 
(7.2-7. S^IO"4  cal/cm-s-K  at  5  atm  up  to  (24-26)- 1  O'4  cal/cm-s-K  at  100  atm.  Values  of  Xs  slightly 
increase  with  T0. 
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Table  3:  Burning  wave  parameters  of  mixture  3  ( [BAMO-AMMOJ/HMX) 


p,  atm 

5 

10 

20 

50 

80 

100 

No 

To,  °C 

20 

-100/20/90 

-100/20/90 

-100/20/90 

-100/20/90 

-100/20 

1 

m, 

g/cm2s 

0.1 

0.16/0.22/ 

0.28 

0.32/0.41/ 

0.48 

0.72/0.85/ 

0.98 

1.06/1.21/ 

1.37 

1.24/ 

1.44 

2 

TS,°C 

330 

345/365/380 

390/407/420 

460/470/480 

483/490/495 

492/500 

3 

cp-10’4,  K/cm 

2.2 

3.4/3.7/5.0 

5.0/6.0/6.5 

8.0/9.0/11 

9.6/11/13.3 

10.4/12 

4 

q,  cal/g 

46 

44.6/37/41 

36/34/33 

27/26/29 

23/24/25 

22/22 

5 

qr,  cal/g 

13.8 

10/8.2/7.7 

7.3/5. 7/6. 2 

3. 8/3. 5/3. 3 

2. 6/2. 5/2.3 

2.4/2. 0 

6 

Q,  cal/g 

83 

142/112/86 

172/136/113 

213/170/143 

227/181/154 

232/187 

7 

/,  nm 

194 

143/111/100 

100/80/77 

46/58/55 

54/50/48 

53/45 

8 

35 

35/35/40 

30/35/40 

25/30/45 

20/30/35 

20/30 

9 

X*  103,  cm2/s 

1.2 

1.3/1. 5/1.7 

1.8/2.0/2.2 

2.8/3. 0/3. 3 

3. 5/3. 7/4. 8 

3. 8/4.0 

10 

ks-iO4, 

cal/cms-K 

7.8 

8.5/9.8/11.0 

11.7/13/14 

18/19.5/21.5 

22.8/24/26 

24.7/26 

11 

Ti,  °C 

900 

-/ 1020/ 1050 

-/-/- 

77- 

77- 

7- 

12 

Tf,°C 

1850 

1950/2000/ 

2100 

2150/2200/ 

2300 

2220/2300/ 

2350 

2220/2320 

2350 

2220/ 

2320 

13 

Li,  urn 

500 

-/400/700 

77- 

77- 

77- 

7- 

14 

L,  mm, 

2.2 

1.7/1. 7/2.0 

1.7/2. 0/2.0 

1.2/1. 8/2.0 

1.1/1. 3/1.5 

1.0/0. 9 

15 

&,  nm 

51 

32/24/20 

17/13/12 

8.3/7. 0/6.2 

5. 8/5. 0/4. 5 

4. 9/4.0 

16 

Xg*l04, 

cal/cm-s-K 

2.1 

2. 1/2. 2/2. 3 

2. 3/2. 3/2.4 

2.5/2. 5/2. 6 

2.6/2.6/2.6 

2.6/2. 6 

17 

0O,  kcal/cm3s 

0.9 

2.2/3.4/5.9 

6.7/10.3/13.0 

24/32.5/45.8 

43/56.6/77.4 

54.8/73 

4.1.6.  Heat  release  in  solid  and  heat  feedback  from  gas  into  solid. 

The  temperature  profiles  and  burning  surface  temperatures  allow  the  important  heat  parameters  of 
the  burning  waves  to  be  estimated.  Heat  flux  from  gas  to  solid  q-m  by  heat  conduction  is  as 
follows: 

q-m=-^g(T)(dT/dx)0; 

Here:  (dT/dx)0=q>  is  the  temperature  gradient  in  gas  close  to  the  burning  surface.  Heat  feedback 
from  gas  into  solid  by  heat  conduction  q  is  as  follows: 

q=-Xg(T)-<p/m; 

Heat  release  in  the  reaction  layer  of  the  condensed  phase  (or  on  the  burning  surface)  Q  is  as  follows: 

Q=c(Ts.-To)-q-qr+qm; 

Where  qr  is  the  radiant  heat  supply  from  flame  into  solid  and  qm  is  heat  of  nitramine  melting. 

It  was  assumed  that  heat  of  HMX’s  melting  is  equal  to  24  cal/g  (per  gramm  of  HMX)  and  that  of 
RDX  is  equal  to  29  cal/g  (per  gramm  of  RDX).  That  implies  that  qm=19.2  cal/g  for  HMX’s 
containing  mixtures  and  qm=23.2  cal/g  for  RDX’s  containing  mixtures. 

The  radiant  heat  supply  was  estimated  by  formula:  qr=k0-kvTf4/m,  where  k0  is  the  coefficient  of 
the  burning  surface  absorption  of  the  flame  radiation  and  -ki  is  the  degree  of  the  flame  blackness.  It 
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was  assumed  that  ko-0.5,  k,«0.1.  It  can  be  seen  that  at  low  pressures  a  significant  radiant  heat 
supply  from  gas  into  solid  is  observed:  qr=(10.4-15)  cal/g  at  5  atm.  However  it  decreases  quickly 
when  pressure  increases  -  at  100  atm  qr  is  equal  to  2-3  cal/g.  Tables  1-4  show  that  values  of  q,  can 
significantly  decrease  at  low  pressures  when  T0  increases. 

Temperature  gradient  cp  in  gas  close  to  the  burning  surface  increases  with  pressure.  It  has  very  large 
values  and  can  be  as  much  as  1.4-105  K/cm.  However  heat  feedback  from  gas  into  solid  by  heat 
conduction  is  relatively  small  because  of  very  low  values  of  XQ.  Tables  1-4  show  that  q  decreases 
when  pressure  increases  -  q  is  equal  to  64-42  cal/g  at  5  atm  and  to  18-26  cal/g  at  100  atm.  Values 

of  q  change  insignificantly  with  T0  variations.  , 

Heat  release  in  solid  Q  increases  with  pressure.  It  has  significant  values  -  65-98  cal/g  at  5  a  m  and 
179-195  cal/g  at  100  atm.  It  is  important  to  stress  that  Q  increases  with  burning  rate,  lables  1-4 
show  that  at  constant  pressures  values  of  Q  significantly  increases  when  T0  decreases.  It  can  be 
explained  by  the  necessaty  of  a  large  heat  consumption  for  the  solid  heating  at  low  sample 

temperatures. 


Table  4:  Burning  wave  parameters  of  mixture  4  ([BAMO-AMMOJ/RDX) 


p,  atm 

5 

10 

20 

50 

80 

100 

No 

oH 

O 

o 

20 

-100/20/90 

-100/20/90 

-100/20/90 

-100/20/90 

-100/20 

1 

m, 

g/cm2s 

0.12 

0.17/0.22/ 

0.33/0.42/ 

0.75/0.87/ 

1.1/1.27/ 

1.47 

1.3/1. 5 

0.27 

0.49 

1.02 

495/500 

2 

r1 

o 

o 

335 

350/370/380 

390/410/420 

460/470/480 

480/490/500 

3 

cp-10-4,  K/cm 

2.5 

3. 6/3. 5/5.0 

6.0/5. 0/7.0 

9.5/8.5/12 

10.5/10.5/14 

11/11 

4 

q,  cal/g 

42 

46.6/36.6/42.6 

41.8/28.6/34.3 

31.6/24.4/30.6 

25/22/25 

22/19 

5 

qr.  cal/g 

10.4 

9/8/8 

6. 5/6/6 

3.4/3.4/3.4 

2.4/2. 5/2.3 

2.0/2. 1 

6 

Q,  cal/g 

98 

148/118/89 

171/145/115 

212/175/145 

228/187/160 

237/194 

7 

/,  nm 

156 

118/107/106 

90/80/74 

63/57/33 

50/47/43 

48/44 

8 

/m,  [im 

90 

60/80/100 

45/65/70 

40/50/60 

40/40/40 

30/44 

9 

2t*103,  cm2/s 

1.2 

1.3/1.4/1.8 

1.9/2.2/2.3 

2.9/3. 1/3.5 

3. 5/3. 9/4. 2 

3. 9/4. 2 

10 

Xs-104, 

7.5 

8.1/8.7/11.2 

11.8/13.7/14 

18/19.3/21.8 

21.8/24/26 

24.3/26 

11 

cal/cms-K 

Ti,°C 

950 

-/1000/1080 

7- 

12 

Tf,  °C 

1800 

1900/2000/ 

2100 

2100/2200/ 

2300 

2200/2300/ 

2400 

2230/2330 

2400 

2230/ 

2330 

13 

Li,  nm 

800 

-/450/600 

-/-/- 

-/-/- 

7- 

14 

L,  mm 

2.5 

1.0/2. 0/2.5 

1. 5/2.0/210 

1.5/1. 5/1. 5 

1.2/1. 3/1.3 

1.0/ 1.0 

15 

&,  urn 

40 

32/25/20 

17/14/12 

8.0/6.8/6.0 

5.6/4.8/4.2 

4.7/4. 1 

16 

XglO4, 

2.0 

2.2/2.3/2.3 

2.3/2.4/2.4 

2.5/2. 5/2.6 

2.6/2.6/2.6 

2.6/2. 6 

17 

cal/cm*sK 
(Po,  kcal/cm3s 

1.2 

2.5/3.2/5.7 

8.3/8.8/14.4 

30/31.4/52 

49/56.7/87.5 

60.8/70 
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4.1.7.  Heat  release  rate  in  gas  close  to  the  burning  surface. 

Heat  release  rates  <t»0  in  gas  close  to  the  surface  was  estimated  by  the  use  of  the  obtained  data.  O0 
shows,  obviously,  the  intensity  of  chemical  reactions  in  gas.  Values  of  O0  were  obtained  from  the 
heat  conduction  equation  and  the  temperature  profiles.  The  heat  conduction  equation  connected 
with  the  burning  surface  for  the  gas  phase  of  stable  propagated  combustion  wave  is  as  follows: 

d/dx-(X.gdT/dx)-m  Cp-dT/dx+0  (x)=0; 

Here  O(x)  is  the  distribution  of  the  heat  release  rate  in  the  gas  phase  of  the  combustion  wave.  The 
first  term  of  the  equation  at  x«0  is  small  and  thus  the  following  approximately  expression  for 
O(x»0)=<Do  is  valid: 

<D0=cp-m<p; 

Tables  1-4  show  that  the  values  00  increase  with  pressure  from  (0.5-0.8)kcal/cm3s  at  2  atm  up  to 
(34-65)  -kcal/cm3s  at  100  atm. 

4.2.  Burning  Rate  Control  Regions. 

The  results  of  measurements  of  q+qr  and  Q  and  the  investigations  of  the  gas  phase  zone 
behaviour  show  that  burning  rate  of  mixtures  1-4  are  caused  by  the  heat  release  in  solid  (or  on  solid 
surface)  and,  in  a  smaller  degree,  by  the  heat  feedback  from  the  gas  layer  which  is  close  to  the 
surface.  It  was  mentioned  above  also  that  the  values  of  of  q+qr  decreases  when  the  burning  rate 
increases.  An  important  conclusion  can  be  made  from  the  analysis  of  the  obtained  data:  the  high 
temperature  region  of  the  burning  waves  cannot  affect  the  burning  rate.  The  influence  of  the  heat 
release  in  the  gas  phase  on  the  burning  rate  can  be  estimated  by  the  following  formula  obtained  as  a 
solution  of  the  heat  conduction  equation: 

oo 

m  q=J  <I>exp(-x/fr)dx;  .  (1.5) 

0 

This  formula  shows  that  influence  of  heat  release  in  gas  on  m  decreases  very  quickly  when  x>$. 
Tables  1-4  show  that  values  &  are  small  and  they  are  significantly  smaller  than  gas  zone  sizes  L  and 
L[.  Indeed,  values  of  fr  are  equal  to  50-4  pm,  but  Li  are  equal  to  200-900  pm  and  L  are  equal  to 
2.5  -1  mm.  The  data  imply  that  high  temperature  region  of  the  burning  waves  cannot  affect  the 
burning  rate. 

Thus  the  obtained  results  show  that  the  burning  rate  control  regions  in  the  combustion  waves  for  all 
regimes  of  the  mixture  combustion  are  the  regions  of  heat  release  in  solid  just  under  the  burning 
surface  (or  immediately  on  the  burning  surface)  and  thin  low-temperature  gas  layers  close  to  the 
burning  surface.  High  temperature  gas  regions  cannot  influence  the  burning  rate  in  fact  because  of  a 
very  large  heat  resistance  of  the  gas  phase  of  the  mixtures. 

4.3.  Macrokinetics  of  Solid  Gasification. 

The  following  equation  connects  burning  rate  of  solid  with  burn-surface  temperature  and 
macrokinetic  characteristics  of  solid  gasification: 

m2=XsP/Q2-9tTs2/Ek0Q*l/N-exp(-E/9tTs); 

where  N=l/ris-l-(l-r)s)/rls-lii(l-r|s)-(q/Q)-Ilii(l-Tis)]/ns;  %=Q/Q*!  Q*  is  the  maximum  of 
the  heat  release  in  solid  E  is  activation  energy  of  limiting  stage  of  the  gasification  process,  Ts  in  K. 
For  studied  combustion  waves  N=l.  The  expression  was  obtained  by  solution  of  system  of  two 
equations  for  the  steady  propagated  burning  wave  -  the  heat  conduction  equation  for  solid  phase 
and  the  equation  of  diffiasion  of  reagents.  The  burning  wave  propagates  due  to  heat  release  in  solid 
Q  and  heat  feedback  from  gas  to  solid  q.  Function  of  the  volumetric  heat  release  rate  Os  in  solid 
was  assumed  as  follows: 

Os(n,  T)=Q*ko  P-(l  -n)  exp(-E/9tT); 
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Here  r|  is  the  reaction  completeness  (n=tis  on  the  burning  surface),  ko  is  the  preexponent 
multiplier. 

The  member  exp(-E/9?Ts)  plays  the  most  important  role  in  this  formula.  Because  of  that  the 
connection  m  and  Ts  can  be  presented  by  the  following  simplified  expression: 

m= Aexp(- E/29tTs) ; 

This  expression  has  been  designated  ‘gasification  law”.  The  results  of  measurement,  presented  in 
Tables  1-4,  allow  the  following  gasification  laws  for  mixtures  1-4  to  be  obtained: 

m=3.8-103-exp(-25000/29lTs);  (m  in  g/cm2s,  T  in  K),  at  TS<470°C 

m=1.7-106-exp(-43000/29lTs);  (m  in  g/cm2s,  T  in  K),  at  TS>470°C 

The  obtained  values  of  activation  energy  of  the  gasification  process  and  the  point  of  E  changing 
(TS*=470°C)  are  very  close  to  those  of  cyclic  nitramines.  Possibly,  gasification  of  oxidizers  (HMX 
and  RDX)  in  the  reaction  layer  of  burning  mixtures  1-4  is  the  limiting  stage  of  the  process  of  the 
mixture  gasification. 

4.4.  Pressure  and  Temperature  Sensitivities  for  m  and  Ts. 

Criteria  of  Stable  Combustion. 

The  following  pressure  and  temperature  sensitivities  for  burning  rate  and  burning  surface 
temperature  were  found  from  the  data  presented  on  Tables  1-4  (T„=20°C): 

P=(51nm/5ToVconst,  -  temperature  sensitivity  of  mass  burning  rate; 
r=(dTJcTTo)p-comt,  -  temperature  sensitivity  of  burning  surface  temperature; 
v=(51nrn/51np)7b-Const,  -  pressure  sensitivity  of  mass  burning  rate; 
p=(Ts-T0)‘l-(5Ts/51np)ro-const,  -  pressure  sensitivity  of  burning  surface  temperature; 

The  obtained  results  are  shown  in  Tables  5  and  6.  Criteria  of  stable  combustion  -  Zeldovich 
criterion  k=P(Ts-T0)  and  Novozilov  criterion  k*=(k- 1  )2/(k+1  )r  -  are  presented  also  in  Table  5  and 
6.  Criterion  k*  is  necessary  to  calculate  only  for  regimes  having  values  of  k>1.  Combustion  is 
stable  if  k*<1 .  Tables  5  and  6  show  that  mixtures  1-4  have  stable  combustion  at  To=20°C. 

At  To=-100°C  exist  regimes  when  k>l,  however  for  these  regimes  always  k*«l. 

At  To=90°C  always  k<l.  Thus  mixtures  1-4  have  stable  combustion  under  investigated  conditions. 


Table  5:  Sensitivities  of  burning  rate  and  surface  temperature 
and  criteria  of  stable  combustion  of  mixtures  1  and  2  at  To=20°C 
(presented  as  fractions:  mix.  l/mix.2) 


p,  atm 

10 

20 

50 

80 

100 

V 

0.83/0.92 

0.80/0.84 

0.76/0.80 

0.76/0.80 

0.76/0.80 

r 

0.17/0.19 

0.155/0.20 

0.13/0.15 

0.10/0.105 

0.09/0.083 

p 

0.172/0.163 

0.171/0.179 

0.115/0.111 

0.104/0.110 

0.092/0.092 

k 

0.97/1.01 

0.87/0.924 

0.79/0.94 

0.79/1.038 

0.8/0.96 

k* 

-/0.0003 

-/- 

-/- 

-/O.OOS 

-/- 
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Table  6:  Sensitivities  of  burning  rate  and  surface  temperature 
and  criteria  of  stable  combustion  of  mixtures  3  and  4  at  To=20°C 
(presented  as  fractions:  mix.3/mix.4) 


p,  atm 

10 

20 

50 

80 

100 

V 

0.90/0.95 

0.82/0.84 

0.80/0.80 

0.79/0.78 

0.79/0.78 

r 

0.19/0.16 

0.16/0.16 

0.12/0.12 

0.066/0.12 

0.067/0.04 

F 

0.174/0.163 

0.168/0.172 

0.098/0.098 

0.094/0.094 

0.092/0.092 

k 

0.95/0.81 

0.76/0.78 

0.72/0.81 

0.61/0.7 

0.576/0.576 

k* 

-A 

7- 

V- 

7- 

7- 

4.5.  Pressure-Driven  Burning-Rate  Response  Functions. 

It  is  interesting  to  investigate  burning  rate  behaviour  of  mixtures  1-4  under  conditions  of  oscillatory 
pressure.  Let  us  assume  that  p=p0+Prcosarf,  where  Pi«p0  and  xu  is  cyclic  frequency  of 
oscillation  in  1/s.  Then  nondimensional  burning  rate  is  equal  to  F=l+Fi-cos crt,  where  V=r\jr^ 
Vi^r^/rbo;  is  a  stationary  burning  rate  and  rh~.  It  is  convenient  to  analyse  burning  rate  response 
in  complex  form.  In  this  case  complex  nondimensional  burning  rate  is  equal  to  V=l+Vve'ia>t  and 
complex  nondimensional  pressure  is  equal  to  77=1+ 771  •e",a*I  where  r|i=Pi/Po,  and  (o  is 
nondimensional  frequency. 

The  following  expression  for  burning  rate  response  function  to  oscillatory  pressure  takes  place: 

U=[v+(v-r-p-k)-(z- 1 )]/[  1  +r-(z-  l)-k-(z- 1  )/z]; 

Here  z=(l+Vl  +  4i<y  )/2;  U=V 1/771,  Vi  is  nondimensional  burning  rate  complex  amplitude,  and  co  is 
nondimensional  frequency  which  is  equal  to  the  cyclic  frequency  w,  in  1/s,  multiplied  by  the  time  of 
solid  heat  layer  relaxation  t>2  Values  of  v,  p,  k  and  r  here  are  based  on  the  mean  pressure  p0  and 
because  of  that  data  from  Tables  5,' 6  can  be  used  for  estimations  of  complex  functions  of  U.  Thus, 
indeed,  we  can  investigate  burning  rate  behaviour  of  mixtures  1-4  under  conditions  of  pulsating 
pressure. 

Results  of  calculationsof  dimensionless  amplitudes  Re{U}  and  phase  Im{U}  of  the  response  as 
functions  of  co  are  shown  on  Figs.  19-22  (T„=20°C).  Response  functions  of  mixtures  1-4  are  close 
each  other  and  because  of  that  Figs.  19-22  show  functions  of  Re{U(<y)}  and  Im{U(<y)}  only  for 
mixture  1.  Figures  19-22  show  that  Re{U}max  are  equal  to  1.3-1. 5  and  practically  independent  on 
pressure.  Resonance  frequencies  co„  are  equal  approximately  to  3.5  at  10  and  20  atm,  to  4.5  at  50 
atm  and  to  6.5  at  100  atm.  It  is  natural  that  function  Im{U(o>)},  presenting  the  phase  of  response 
functions,  changes  the  sign  at  the  resonance  frequency.  Functions  of  Re{U((y)}  weakly  depend  on 
0),  especially  at  elevated  p.  The  result  shows  that  resonance  phenomena  in  burning  mixtures  1-4  are 
significantly  depressed  due  to  viscous  behaviour  of  the  oscillations.  In  fact,  the  pulsating 
combustion  of  the  mixtures  is  an  example  of  a  relatively  smoothed  resonance  phenomenon.  Such 
behaviour  must  be  attributed  to  large  amount  of  cyclic  nitramines  in  the  mixtures,  because  these 
nitramines  have  decreased  Re{U}max  and  weak  dependence  Re{  U(o») } . 
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5.  Results  of  Measurements  and  Processing  for  Mixtures  5-8. 

The  characteristics  of  the  initial  mixtures  were  as  follows: 

5)  GAP1U/HMX,  density  p=l,38  g/cm3;  sum  formula:  O23.6H31.7N27.6C16.9;  Tad=2695  K; 

6)  GAP1U/RDX,  density  p=  1,19  g/cm3;  sum  formula:  O23.6H31.7N27.66Ci6.86;  Tad=2705  K; 

7) (BAMO-THF)/HMX,  density  p=l,70  g/cm3;  sum  formula:  O23.3H34.9N26.6C18.3;  Tad=2465  K; 

8)  (BAMO-THF)/RDX,  density  p=l,65  g/cm3;  sum  formula:  O23.4H33.9N26.1C18.7;  Tad=2445  K; 
Binder  (BAMO-THF)  have  proportion  50:50.  Cylindrical  samples  of  8  mm  in  diameter  and  2  cm  in 
length  were  burned  in  a  bomb  of  constant  pressure.  Temperature  profiles  and  surface  temperatures 
were  obtained  for  the  indicated  mixtures  at  pressures  10-100  atm  and  at  sample  temperature 
To=+20C°.  Carbon  residue  on  the  burning  surface  of  the  mixtures  did  not  observed.  Temperature 
profiles  were  received  and  after  averaging  they  were  analysed.  Tables  7,  8  show  the  obtained 
characteristics.  The  nominations  are  the  same  as  in  Tables  1-4. 

5. 1 .  Burning  Wave  Characteristics. 

5.1.1.  Mass  burning  rate.  Mass  burning  rate  m  increases  for  mixtures  5-8  when  pressure  increases. 
Values  of  m  for  mixtures  5-7  are  close  each  other  and  that  for  mixture  8  is  less,  especially  at 
elevated  pressures.  At  2  atm  m  are  equal  to  0.06-0.08  g/cm2s,  at  10  atm  -  0.17-0.22  g/cm2s,  at  20 
atm  -  0.3-0.45  g/cm2s.  At  elevated  pressures  diapason  of  values  m  increases:  at  50  atm  -  0.8-0.82 
g/cm2s  for  mixtures  6  and  7,  but  for  mixture  8  -  0.45  g/cm2s;  at  100  atm  -  1.38-1.65  g/cm2s  for 
mixtures  6  and  7,  but  for  mixture  8-1.09  g/cm2s. 

5.1.2.  Burning  surface  temperature.  Burning  surface  temperature  Ts  always  increases  when  m  and 
pressure  increase.  Values  of  Ts  changes  significantly:  from  280-360°C  at  2  atm  up  to  560-580°C  at 
100  atm. 

5.1.3.  Geometric  structure  of  the  temperature  profiles  in  gas.  The  gas  phase  has  two-zone  structure: 
the  first  zone  with  mean  zone  temperatures  Ti  and  the  flame  zone  with  the  final  zone  temperatures 
Tf.  Both  zones  merge  at  20-50  atm.  Only  the  first  zone  exists  at  2  atm  for  mixture  6  with  the  final 
zone  temperature  Ti.  Values  of  Ti  are  equal  to  900-1 100°C.  Thermodynamically  calculated  values 
of  Tf  are  achieved  at  20-50  atm.  Distances  Li  between  the  burning  surface  and  the  beginning  of  the 
flame  zone  at  pressures  2  (10)-20  atm  are  very  small  for  mixtures  5,  6  (250-100  pm)  and  very  large 
for  mixtures  7,  8  (0.7-1. 5  mm).  Distances  L  between  the  burning  surface  and  the  flame  -  1.5  -  0.8 
mm  for  mixtures  5,  6  and  3.5  -  1.0  mm  for  mixtures7,  8. 

5.1.4.  Characteristic  sizes  in  gas  and  solid.  Mixtures  5,  6  have  /  =210-30  pm  and  mixtures  7,  8  -  / 
=65-20  pm.  Conductive  sizes  $  of  the  gas  phase  decrease  when  pressure  increases  and  comprise 
23-3  pm  for  mixtures  7,  8  and  85-4.8  pm  for  mixtures  5,  6.  It  is  obvious  that  ft«Li. 

5.1.5.  Heat  diffusivity  and  heat  conductivity  of  heat  layer  in  solid.  Values  x  increase  from  (1- 
1.5)- 1 0'3  cm2/s  at  low  pressures  up  to  (3-3.5)- 10'3  cm2/s  at  100  atm.  That  are  normal  values  of  %  for 
many  propellant  mixtures.  The  solid  heat  conductivity  was  obtained  by  formula  A.s=x-c  p;  it  was 
assumed  that  solid  specific  heat  c=0.4  cal/g-K.  Tables  7,  8  show  that  h  increases,  when  pressure 
increases,  from  510'4cal/cm-s  K  at  2  atm  up  to  (16-23)10'4cal/cm-s  K  at  100  atm. 

5.1.6.  Heat  release  in  solid  and  heat  feedback  from  gas  to  solid.  The  radiant  heat  supply  was 
estimated  by  formula:  qr=ko-ki-Tf4/m,  where  k0  is  coefficient  of  the  burning  surface  absorption  of 
the  radiation  and  -ki  is  degree  of  the  flame  blackness.  That  was  assumed  that  ko=0.5,  ki=0.1.  It  can 
be  seen  that  at  low  pressures  a  significant  radiant  heat  supply  from  gas  into  solid  is  observed. 
Tables  7,  8  show  that  the  temperature  gradient  cp  in  gas  close  to  the  burning  surface  has  very  large 
values:  it  can  be  as  much  as  10s  K/cm.  Values  of  q  and  Q  are  comparable  at  2-5  atm.  But  at 
elevated  pressures  values  Q  are  more  ,  than  values  of  q.  As  a  rule,  Q  increases  with  pressure, 
however  q,  as  a  rule,  decreases. 
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5.1.7.  Heat  release  rate  in  gas  close  to  the  burning  surface.  Tables  7,  8  show  that  values  of  <D, 
increase  with  pressure  from  (0.5-0.8)-kcal/cm3s  at  2  atm  up  to  (34-65)  •kcal/cm3s  at  100  atm. 


Table  7:  Burning  wave  parameters  of  mixtures  5  and  6  at  To=20°C 
(indicated  as  fractions:  mix.5/mix.6) 


p,  atm 

2 

5 

10 

20 

50 

100 

fbi  cm/s 

0.06/0.047 

0.107/0.126 

0.17/0.22 

0.30/0.38 

0.58/- 

1.0/- 

m,  g/cm2s 

0.083/0.056 

0.148/0.15 

0.235/0.27 

0.42/0.455 

0.8/- 

1.38/- 

r* 

o 

O 

360/280 

400/370 

440/440 

490/460 

560/- 

580/- 

E 

o 

5? 

o 

9- 

2. 3/2.0 

4.5/3. 5 

5.0/4. 5 

6.5/6.0 

10/- 

ii/- 

q,  cal/g 

55/68 

•  70/49 

51/40 

39/32 

34/- 

22/- 

Q,  cal/g 

74/41 

84/100 

124/142 

162/161 

198/- 

219/- 

/,  nm 

170/210 

110/80 

90/70 

60/45 

35/- 

30/- 

Tl,  °C 

950/900 

1000/1000 

1050/1100 

1100/- 

-/- 

7- 

Tf,°C 

2100/1700 

2200/2100 

2250/2200 

2300/2300 

2400/- 

2420/- 

Li,  Jim 

250/- 

150/230 

120/180 

100/- 

-/- 

7- 

L,  mm 

1.5/1. 5 

1.0/1. 2 

1.0/ 1.0 

1.2/0. 9 

0.8/- 

0.8/- 

S,  urn 

63/85 

84/35 

24/21 

14/13 

8/- 

4.8/- 

X*103,  cm2/s 

1/1 

1.2/1. 2 

1.5/1. 5 

1. 8/2.0 

2/- 

3/- 

Imi 

60/100 

50/60 

40/55 

35/40 

30/- 

25/- 

qr,  cal/g 

26/18 

17/14 

12/9 

7/6.5 

4/- 

2/- 

A,s-104,  cal/cm  s-K 

5.5/4. 8 

'  -6. 6/5.8 

8. 3/7.2 

9.9/9.6 

11/- 

16.5/- 

X,g-104,  cal/cm  s-K 

2/1.9 

2.3/2. 1 

2A/2A 

2. 5/2.4 

2.7/- 

2.8/- 

Oo,  kcal/cm3s 

0.8/0.47 

2.8/2.2 

4.9/5. 1 

11.5/11.5 

34.0/- 

64.5/- 

5.2.  Burning  Rate  Control  Regions. 

The  results  of  measurements  of  q+qr  and  Q  and  the  investigations  of  the  gas  phase  zone 
behaviour  show  that  burning  rate  of  mixtures  5-8  are  caused  by  the  heat  release  in  solid  (or  on  solid 
surface)  and,  in  a  smaller  degree,  by  the  heat  feedback  from  the  gas  layer  which  is  close  to  the 
surface.  An  important  conclusion  can  be  made  from  the  analysis  of  the  obtained  data:  the  high 
temperature  region  of  the  burning  waves  cannot  affect  the  burning  rate.  The  same  analysis  as  in  the 
previous  case  show  that  high  temperature  region  of  the  burning  waves  cannot  affect  the  burning 
rate.  Thus  the  obtained  results  show  that  the  burning  rate  control  regions  in  the  combustion  waves 
for  all  regimes  of  the  mixture  combustion  are  the  regions  of  heat  release  in  solid  just  under  the 
burning  surface  (or  immediately  on  the  burning  surface)  and  thin  low-temperature  gas  layers  close 
to  the  burning  surface.  High  temperature  gas  regions  cannot  influence  the  burning  rate  in  fact 
because  of  a  very  large  heat  resistance  of  the  gas  phase  of  the  mixtures. 
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Table  8:  Burning  wave  parameters  of  mixtures  7  and  8  at  To=20°C 
(indicated  as  fractions:  mix.  7/mix.  8) 


p,  atm 

10 

20 

50 

100 

rb(  cm/s 

0.11/0.14 

0.18/0.2 

0.26/0.5 

0.64/1.0 

m,  g/cm2s 

0.187/0.23 

0.306/0.33 

0.45/82 

1.09/1.65 

r1 

o 

O 

380/400 

420/450 

450/470 

500/520 

(p-KH,  K/cm 

2. 8/2.5 

3.0/4.0 

6. 5/7. 3 

7. 5/7. 5 

q,  cal/g 

27/19.6 

17.6/23 

27.4/17 

13.8/9 

Q,  cal/g 

128/148 

156/166 

158/183 

195/213 

/,  urn 

64/65 

65/70 

56/37 

32/20 

Ti,°C 

950/1000 

1000/1100 

-/- 

-/- 

Tf,  °C 

1880/2000 

1950/2100 

2200/2170 

2200/2170 

Li,  mm 

1. 5/1.0 

0.8/0.7 

7- 

7- 

L,  mm 

3. 5/2.5 

2. 5/1. 5 

2.0/2. 0 

1.0/1. 5 

8,  pm 

23/19 

14/14 

10/5.5 

4.3/3.0 

X*103,  cm2/s 

1.2/1. 5 

2.0/2. 3 

2.5/3.0 

3.5/3. 5 

/m,  jam 

30/50 

55/60 

40/30 

20/20 

qr,  cal/g 

7. 8/7.9 

5. 5/6. 5 

5. 6/2.9 

2.3/1. 5 

A,s*104,  cal/cm-s-K 

8. 2/9. 5 

13.6/15.2 

17/19.8 

23.8/23.0 

Xg-104,  cal/cm-s-K 

1.8/1. 8 

1.8/1. 9 

1.9/1. 9 

2.0/2.0 

Oo,  kcal/cm3s 

2. 2/2.4 . 

3.8/5. 5 

12.3/26.2 

34.6/52.6 

6.  Conclusions  About  Combustion  Mechanisms  of  Nitramine/Polymer  Mixtures. 

Nitramine/Energetic  Polymer  Binders. 

1)  The  mixtures  1-8  have  large  heat  release  in  solid  Q,  which  increases  with  pressure  and  burning 
rate  m.  Heat  feedback  from  gas  to  solid  is  less  (than  Q)  and  decreases  when  p  and  m  increase. 

At  low  pressures  a  significant  role  in  the  heat  feedback  from  gas  to  solid  plays  radiant  heat  flux 
from  flame  to  solid. 

2)  Burning  rate  control  regions  are  located  on  solid  surface  (due  to  large  heat  release  in  solid)  and 
in  thin  low-temperature  gas  layer  close  to  the  burning  surface. 

3)  Process  of  oxidizer  decomposition  plays  the  limiting  role  in  the  mixtures  gasification,  because 
macrokinetics  of  the  mixture  gasification  are  very  close  to  those  of  HMX  and  RDX. 

4)  The  gas  phase  of  burning  mixtures  has  two-zone  structure  at  low  pressures  (5-10  atm),  which 
merges  into  one  zone  when  pressure  increases.  Increasing  initial  temperature  can  create  two-zone 


14 


structure  from  one  zone  and  decreasing  initial  temperature  can  merge  the  two-zone  structure  into 
one  zone. 

5)  Thickness  L  of  the  reacting  gas  zone  is  equal  to  1-2.5  mm.  As  a  rule,  L  decreases  when  pressure 
increases  (especially  in  region  50-100  atm)  and  increases  with  the  initial  temperature. 

6)  Completeness  of  combustion  is  achieved  at  50-80  atm.  Flame  temperature  increases  with  the 
initial  temperature. 

7)  Stable  combustion  is  observed  under  investigated  conditions. 

8)  Pressure-driven  bum-rate  response  functions  of  the  investigated  mistures  have  viscous 
behaviour,  without  charp  increase  of  burning  rate  response  amplitude  close  to  resonance  frequency. 

Comparison  of  Combustion  Mechanisms  of  Mixtures  based  on  Energetic  Binders  and 

Mixtures  based  on  HTPB. 

1)  Mixtures  based  on  HTPB  have  decreased  burning  rate,  burning  surface  temperature,  heat  release 
in  solid  and  heat  feedback  from  gas  to  solid,  than,  under  comparible  conditions  (pressure,  initial 
temperature),  mixtures  based  on  energetic  binders. 

2)  Mixtures  based  on  HTPB  have  increased  thickness  of  reacting  gas  zone.  Two-zone  structures  in 
the  gas  phase  are  observed  up  to  higher  pressures  than  two-zone  structures  of  mixtures  based  on 
energetic  binders. 

3)  Mixtures  based  on  HTPB  have  decreased  flame  temperatures. 

4)  Mixtures  based  on  HTPB  and  on  energetic  binders  have  similar  macrokinetics  of  solid 
gasification,  and  similar  positions  of  burning  rate  control  regions. 
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FIGURES 


-0.5  0  0.5  1  1.5  2 

distance  x,  mm 

Figure  1:  Averaged  temperature  profiles  T(x)  of  mixtures  1 
(GAP2/HMX)  and  2  (GAP2/RDX)  at  5  atm,  To=20'C. 


-0.5  0  0.5  1  1.5  2  2.5 

distance  x,  mm 

Figure  2:  Averaged  temperature  profiles  T(x)  of  mixtures  3 
([BAM0-AMM0]/HMX)  and  4  ([BAMO-AMMO]/RDX) 
at  5  atm,  To=20’C. 
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temperature  T’C  temperature  T'C 


Figure  3:  Averaged  temperature  profiles  T(x)  of  mixture  1 
(GAP2/HMX)  at  10  atm  and  different  To. 


Figure  4:  Averaged  temperature  profiles  T(x)  of  mixture  1 
(GAP2/HMX)  at  20  atm  and  different  To. 
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temperature  T'C  temperature  T’C 


Figure  5:  Averaged  temperature  profiles  T(x)  of  mixture  1 
(GAP2/HMX)  at  50  atm  and  different  To. 


Figure  6:  Averaged  temperature  profiles  T(x)  of  mixture 
1  (GAP2/HMX)  at  100  atm  and  different  To. 
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temperature  T’C  temperature  T'C 


Figure  7:  Averaged  temperature  profiles  T(x)  of  mixture  2 
(GAP2/RDX)  at  10  atm  and  different  To. 


Figure  8:  Averaged  temperature  profiles  T(x)  of  mixture  2 
(GAP2/RDX)  at  20  atm  and  different  To. 
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Figure  9:  Averaged  temperature  profiles  T(x)  of  mixture  2 
(GAP2/RDX)  at  50  atm  and  different  To. 


Figure  10:  Averaged  temperature  profiles  T(x)  of  mixture  2 
(GAP2/RDX)  at  100  atm  and  different  To. 
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Figure  11:  Averaged  temperature  profiles  T(x)  of  mixture  3 
[(BAMO-AMMO)/HMX]  at  10  atm  and  different  To. 


Figure  12:  Averaged  temperature  profiles  T(x)  of  mixture  3 
[(BAM0-AMM0)/HMX]  at  20  atm  and  different  To. 
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Figure  13:  Averaged  temperature  profiles  T(x)  of  mixture  3 
[(BAMO-AMMO)/HMX]  at  50  atm  and  different  To. 


Figure  14:  Averaged  temperature  profiles  T(x)  of  mixture  3 
[(BAMO-AMMO)/HMX]  at  100  atm  and  different  To. 
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Figure  15:  Averaged  temperature  profiles  T(x)  of  mixture  4 
[(BAMO-AMMO)/RDX]  at  10  atm  and  different  To. 
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Figure  16:  Averaged  temperature  profiles  T(x)  of  mixture  4 
[(BAM0-AMM0)/RDX]  at  20  atm  and  different  To. 
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Figure  17:  Averaged  temperature  profiles  T(x)  of  mixture  4 
[(BAMO-AMMO)/RDX]  at  50  atm  and  different  To. 


Figure  18:  Averaged  temperature  profiles  T(x)  of  mixture  4 
[(BAMO-AMMO)/RDX]  at  100  atm  and  different  To. 
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